February 2, 2008 5:11 Proceedings Trim Size: 9in x 6in 



schien'ring05 



O 
O 

(N ; HEAVY-FLAVOR PHOTO- AND ELECTROPRODUCTION 

X> 1 AT NLO 

r-i ■ 

I. SCHIENBEIN 

Southern Methodist University 
C*") ' Department of Physics 

^ , 104 Fondren Science Building 

3215 Daniel Avenue 
Dallas, Texas 75275-0175, USA 
E-mail: schien@physics.smu.edu 



O 

o 



We review one-particle inclusive production of heavy-flavored hadrons in a frame- 
work which resums the large collinear logarithms through the evolution of the FFs 
^ (H and PDFs and, at the same time, retains the full dependence on the heavy-quark 

O ' mass without additional theoretical assumptions. The main focus is on the pro- 

I ' duction of D* mesons in deep inelastic electron— proton scattering at HERA. We 

show results, neglecting for the time being the heavy-quark mass terms, for deep 
inelastic D* meson production at finite transverse momenta. Work to implement 
i"^H ' this process in the above mentioned massive QCD framework is in progress. 
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1. Introduction 

One-particle inclusive production processes provide extensive tests of per- 
turbative quantum chromodynamics (pQCD). In contrast to fully inclusive 
processes, it is possible to study distributions in the momentum of the final 
state particle and to apply kinematical cuts which come close to the ex- 
perimental situation. On the other hand, contrary to even more exclusive 
cases, QCD factorization theorems 1,2 still hold stating that this class of ob- 
servables can be computed as convolutions of universal parton distribution 
functions (PDFs) and fragmentation functions (FFs) with perturbatively 
calculable hard scattering cross sections. As is well-known, it is due to the 
factorization property that the parton model of QCD has predictive power. 
Hence, tests of the universality of the PDFs and FFs are of crucial impor- 
tance for validating this QCD framework. At the same time, lowest order 
expressions for the hard scattering cross sections are often not sufficient for 
meaningful tests and the use of higher order computations is mandatory. 
The perturbative analysis is becoming more involved and interesting 
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if the observed final state hadron contains a heavy (charm or bottom) 
quark. In this case, the heavy-quark mass mh enters as an additional 
scale. Clearly, the conventional massless formalism, also known as zero- 
mass variable-flavor- number scheme (ZM-VFNS), can also be applied to 
this case, provided the hard scale Q of the process is much bigger than 
the heavy-quark mass such that terms nih/Q are negligible. However, at 
present collider energies, most of the experimental data lie in the kinematic 
region Q > mh and it is necessary to take the power-like mass terms into 
account in a consistent framework. 

The subject of this article is the theoretical description of one-particle 
inclusive production of heavy-flavored hadrons Xh = D,B,A C ,... in a 
massive variable-flavor- number scheme (GM-VFNS). In such a scheme the 
large collinear logarithms of the heavy-quark mass ln/i/m/j are subtracted 
from the hard scattering cross sections and resummed through the evolu- 
tion of the fragmentation functions (FFs) and parton distribution functions 
(PDFs). At the same time, finite non-logarithmic mass terms rrih/Q are 
kept in the hard part and fully taken into account. 

In order to test the pQCD formalism, in particular the universality of 
the FFs, it is important to provide a description of all relevant processes 
in a coherent framework. Therefore, it is important to work out the GM- 
VFNS at next-to-leading order (NLO) of QCD for all the relevant processes. 
Previously, the GM-VFNS has been applied to the following processes: 
7 + 7 — ► D* + + X (direct part) 3 , 7 + 7 — > D* + + X (single resolved part) 
4 , 7 +p -» D*+ + X (direct part) 5 , p + p -» (L>° , D*+ , D+ , D+) + X 
6 - 7 ' 8 , where the latter results for hadron-hadron collisions also constitute 
the resolved contribution to the photoproduction process 7 + p — > Xh + X. 

In this contribution, we will review the production of heavy-flavored 
hadrons Xh in hadron-hadron, photon-proton and deep inelastic electron- 
proton collisions where the main focus will be on the electroproduction case. 
We will show results, neglecting for the time being the rrih/Q mass terms, 
for deep inelastic D* meson production at finite transverse momenta. 

2. Theoretical Framework 
2.1. GM-VFNS 

The differential cross sections for inclusive heavy-flavored hadron produc- 
tion can be computed in the GM-VFNS according to the familiar factor- 
ization formulae, however, with heavy-quark mass terms included in the 
hard scattering cross sections 9 . Generically, the physical cross sections are 
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expressed as convolutions of PDFs for the incoming hadron(s), hard scat- 
tering cross sections, and FFs for the fragmentation of the outgoing parton 
into the observed hadron. All possible partonic subprocesses are taken into 
account. The massive hard scattering cross sections are constructed in a 
way that in the limit m,h — > the conventional ZM-VFNS is recovered. 
A more detailed discussion of the GM-VFNS and the construction of the 
massive hard scattering cross sections can be found in Refs. 6, 7 and the 
conference proceedings 10 > n . 

2.2. Fragmentation Functions 

A crucial ingredient entering these calculation are the non-perturbative FFs 
for the transition of the final state parton into the observed hadron Xh- For 
charm-flavored mesons, X c , three sets of FFs have been employed in our 
analyses: (i) BKK98-D 12 : For X c = D*+ , such FFs were extracted at LO 
and NLO in the MS factorization scheme with n/ = 5 masslcss quark flavors 
several years ago 12 from the scaled-energy (x) distribution da/dx of the 
cross section of e+e~ — > D* + + X measured by the ALEPH 13 and OPAL 
14 Collaborations at CERN LEP1. (ii) KK05-D 15 : Recently, the BKK98- 
D analysis was extended 15 to include X c = D°, D + , £)+, A+ by exploiting 
appropriate OPAL data 16 . (iii) KK05-D2: In Refs. 12 ' 15 , the starting scales 
Uq for the DGLAP evolution of the a — > X c FFs in the factorization scale u' F 
have been taken to be uo = 2m c for a = g, u,u, d, d, s,s, c,c and uq = 2nib, 
with m b — 5 GeV, for a = b,b. The FFs for a = g,u,u,d,d,s,s were 
assumed to be zero at a' F = ^ and were generated through the DGLAP 
evolution to larger values of u' F . For consistency with the MS prescription 
for PDFs, we repeated the fits of the X c FFs for the choice fi — m c ^6- 
This changes the c-quark FFs only marginally, but has an appreciable effect 
on the gluon FF, which is important at Tevatron energies, as was found for 
D* + production in Ref. 6 . These new FFs will be presented elsewhere. 

2.3. Input Parameters 

If not stated otherwise, the following parameters have been chosen for the 
numerical results presented below. For the proton PDFs we have employed 
the CTEQ6M/CTEQ6.1M PDFs from the CTEQ Collaboration 17 and the 
KK05-D2 FFs. We have set m c = 1.5 GeV, m b = 5 GeV and have used the 
two-loop formula for ai™^ (/xr) in the MS scheme with af (mz) = 0.118. 
Our theoretical predictions depend on three scales, the rcnormalization 
scale Hr, and the initial- and final-state factorization scales u F and u' F , 
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Figure 1. Comparison of the CDF data 18 with our NLO predictions for _D*+ . The 
solid line represents our default prediction obtained with /i^ = /ip = /i' F = tut, while 
the dashed lines indicate the scale uncertainty estimated by varying fiR, fip, and jx' F 
independently within a factor of 2 up and down relative to the central values. The right 
figure shows the data-over-theory representation with respect to our default prediction. 

respectively. Our default choice for hadro- and photoproduction has been 
= jiF = \i!p = tut, where ttit = \/~Vt + ™^ is the transverse mass. The 
scale choice in the electroproduction case will be specified below. 

3. Hadroproduction 

Recently, the CDF collaboration has published first cross section data for 
inclusive production of D°, D + , and Df mesons in pp collisions 18 

obtained in Run II at the Tevatron at center-of-mass energies of VS= 1.96 
TeV. The data come as distributions da/dpx with y integrated over the 
range \y\ < 1 and the particle and antiparticlc contributions are averaged. 

Our theoretical predictions in the GM-VFNS are compared with the 
CDF data for D* mesons on an absolute scale in Fig. 1 (left) and in the 
data-over-theory representation with respect to our default results in Fig. 1 
(right). We find good agreement in the sense that the theoretical and 
experimental errors overlap where the experimental results are gathered on 
the upper side of the theoretical error band, corresponding to a small value 
of p,R and large values of hf and [i' F , the [ir dependence being dominant 
in the upper px range. As is evident from Fig. I (right), the central data 
points tend to overshoot the central QCD prediction by a factor of about 1.5 
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at the lower end of the considered pr range, where the errors are largest, 
however. This factor is rapidly approaching unity as the value of pr is 
increased. The tendency of measurements of inclusive hadroproduction in 
Tevatron run II to prefer smaller renormalization scales is familiar from 
single jets, which actually favor = pt/2 19 ■ 

For more details and a comparison with the data for the D°, D + , and 
Df mesons we refer to Ref. 8. A comparison of NLO predictions for p+p — > 
B+ + X in the GM-VFNS with recent CDF data 20 is in preparation 21 . 

4. Photoproduction 

Inclusive photoproduction of D* mesons, 7 + p — > D* + X, has been stud- 
ied in Ref. 5 where the direct part has been computed in the GM-VFNS 
whereas the resolved part has been included in the ZM-VFNS. In this analy- 
sis the BKK98-D FFs have been utilized and, for the resolved contribution, 
the GRV92 photon PDFs 22 . The other parameters have been chosen as 
specified in Sec. 2.3. In Fig. 6 of Ref. 5, the central numerical predictions 
for the transverse momentum (pt) distributions of the D* meson have been 
compared with preliminary ZEUS data 23 . There exist similar data by the 
HI collaboration 24 which have not been used in this analysis. As can be 
seen in this figure, the agreement of the p-r-distributions with the data is 
quite good down to px — 2m c and the mass effects turn out to be small. 
In order to extend the range of applicability of the GM-VFNS into the re- 
gion pt < 3 GeV more work on the matching to the 3-fixcd flavor theory 
would be needed. The Figs. 7 - 9 of Ref. 5, showing results for the rapidity 
(y), invariant mass (W) and inelasticity (z(D*)) distributions, have to be 
taken with a grain of salt since they receive large contributions from the 
transverse momentum region 1.9 < pr < 3 GeV which is outside the range 
of validity of the present theory. 

With the work in Ref. 6, it is now possible to include also the resolved 
part in the GM-VFNS. It will be interesting to compare the complete GM- 
VFNS framework at NLO of QCD, combined with updated FFs, in more 
detail with ZEUS and HI photoproduction data once they are finalized. 

5. Electroproduction 

Recently, the single inclusive electroproduction of light hadrons at finite 
transverse momenta has attracted quite a lot of interest, where the outgoing 
hadron is required to carry a non- vanishing transverse momentum (p T ) in 
the center-of-mass system (CMS) of the virtual photon and the incoming 
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proton. The following partonic subprocesses contribute at leading order 
(LO) and are of the order 0(a s ): 7* + q — > q + g and 7* + g — > g + 
(j. Very recently, the NLO (0(a 2 )) corrections to this process have been 
accomplished by three independent groups 25 > 26 > 27 . Suffice to say here that 
the NLO corrections increase the LO results in certain kincmatical regions 
by large factors and are essential for bringing theory in agreement with the 
experimental results. For more details see Ref. 28. 

Endowed with appropriate fragmentation functions, the computation 
in Ref. 27 has been employed to obtain predictions in the ZM-VFNS for 
the production of heavy-flavored hadrons in electron-proton collisions at 
HERA, e+p — > e + Xh + X. The electron and proton energies have been set 
to E e = 27.5 GeV and E p = 820 GeV (Fig. 2) resp. 920 GeV (Figs. 3, 4) in 
the laboratory frame. Furthermore, the renormalization and factorization 
scales have been chosen 27 as p? R = fip = n' F 2 = .e^L+ipiL xhc dimen- 
sionless parameter £ has been varied between 1/2 and 2 about the default 
value 1 in order to estimate the scale uncertainties of the theoretical pre- 
dictions. The other input parameters have been chosen as described in Sec. 
2.3. The following cuts have been imposed on the numerical results in Figs. 
2-4: 2 < Q 2 < 100 GeV 2 , 0.05 < y < 0.7, 1.5 < p T ^b{D*) < 15 GeV, 
and |??Lab(-D*)| < 1.5 where the momentum transfer Q 2 , the inelasticity y, 
and the pscudorapidity 77 are defined as usual. Furthermore, in Figs. 3-4 
we have asked for the additional constraint p^{D*) > 2 GeV where p^(D*) 
is the transverse momentum of the D* meson in the 7*p-CMS. This cut is 
essential to avoid the collinear singularities as p* T — ► 0. The results have 
been calculated with nj = 5 flavors and include the contribution where a 
bottom quark fragments into the D* meson via b — > B — > D*. No dis- 
tinction is made between D* + and D*~ . Figure 2 shows the p* T spectrum 
in comparison with HI data 29 , collected in the years 1994 to 1996 with 
27.5 GeV positrons colliding with 820 GeV protons at CMS energies of 
yfS = 300 GeV. The data points are reasonably well described if one keeps 
in mind that the theory is expected to work for larger p^ > 2 GeV. Fur- 
thermore, in the ZM-VFNS the heavy-quark mass terms are missing which 
are potentially important in the region of small p* T . These terms will be in- 
cluded as soon as ongoing work to implement this process in the GM-VFNS 
has been completed. 

Finally, results for other distributions (pr.Lab, ?7Lab, W, z, Q 2 , and xsj) 
are presented in Figs. 3-4. 
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Figure 2. NLO QCD predictions in the ZM-VFNS for electroproduction of D* mesons 
in dependence of the transverse momentum of the D* meson in the 7*p CMS. HI data 
29 are shown for comparison. See the text for further details. 
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Figure 3. NLO QCD predictions in the ZM-VFNS for electroproduction of D* mesons 
in dependence of the transverse momentum (left) and the pseudo-rapidity (right) of the 
D* meson in the laboratory frame. For further details on the kinematical cuts and the 
uncertainty band see the text. 



6. Summary 

We have discussed one-particle inclusive production of heavy-flavored 
hadrons in hadron-hadron, photon-proton, and electron-proton collisions 
in a massive variable-flavor- number scheme (GM-VFNS). The importance 
of a unified treatment of all these processes, based on QCD factorization 
theorems, has been emphasized, in order to provide meaningful tests of the 
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Figure 4. As in Fig. 3 for the 7*p invariant mass W, the inelasticity z, the momentum 
transfer Q 2 and Bjorken-x. 



universality of the FFs and hence of QCD. At the same time, it is necessary 
to incorporate heavy-quark mass effects in the formalism since many of the 
present experimental data points lie in a kinematical region where the hard 
scale of the process is not much larger than the heavy-quark mass. This 
is achieved in the GM-VFNS, which includes heavy-quark mass effects in 
a rigorous way and still relies on QCD factorization. We have discussed 
numerical results for the three reactions. In general, the description of the 
transverse momentum spectra is quite good down to transverse momenta 
Pt — 1m%. Extending the range of applicability of our scheme to smaller px 
would require more work on the matching to the corresponding theories in 
the fixed flavor number schemes. Our ZM-VFNS results for the electropro- 
duction of D* mesons indicate that future experimental results by the HI 
collaboration can be nicely described. All this leads us to the expectation 
that a good overall description of the data can be reached in the future 
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